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D Abstract 

In the Next-to-Minimal Supersymmetric Standard Model, a bino-like LSP can be 
as light as a few GeV and satisfy WM AP constraints on the dark matter relic density 
in the presence of a light CP-odd Higgs scalar. We study upper bounds on the direct 
detection cross sections for such a light LSP in the mass range 2 — 20 GeV in the 
NMSSM, respecting all constraints from B-physics and LEP. The OPAL constraints 
on e + e~ — > XiX? (* > 1) play an important role and are discussed in some detail. 
C ■ ' The resulting upper bounds on the spin-independent and spin-dependent nucleon 



cross sections are ~ 10 cm and ~ 4 x 10 cm" , respectively. Hence the 
upper bound on the spin-independent cross section is below the DAMA and CoGeNT 
regions, but could be compatible with the two events observed by CDMS-II. 



1 Introduction 

The DAMA [1] and CoGeNT [2] dark matter detection experiments have reported events in 
excess of the expected background, which would be compatible with a WIMP mass of a few 
GeV. Also the CDMS-II experiment [3] has reported two events, which could be explained 
by a WIMP mass of > 10 GeV (or background). On the other hand, exclusion limits 
from the XenonlO [4], XenonlOO [5] and CDMS-Si [6] experiments set upper bounds on 
the spin independent detection cross sections for this mass range of a WIMP, which seem 
incompatible with the reported hints for a signal. 

In any case, it is important to know whether specific models for dark matter with a 
WIMP mass of a few GeV can produce direct detection cross sections compatible with the 
reported excesses, and/or whether regions in the parameter space of such models can be 
tested by present and future exclusion limits. 

Supersymmetric (Susy) extensions of the Standard Model are popular, amongst others, 
since they predict naturally (for unbroken R-parity, and for a neutral Lightest Supersym- 
metric Particle (LSP)) a candidate for dark matter, with a relic density compatible with 
WMAP constraints [7]. Within the Minimal Supersymmetric extension of the Standard 
Model (MSSM) four neutral fermions (neutralinos x°n i = 1---4) exist, which are com- 
posed of the bino (superpartner of the U(l)y gauge boson), the wino (superpartner of the 
W^ gauge boson) and two higgsinos (superpartners of neutral Higgs bosons). These states 
mix, and the lightest neutralino Xi, which is the lightest eigenvalue of the 4x4 mass matrix, 
will be the LSP (leaving aside the possibility of a sneutrino LSP). 

Often the LSP is dominantly bino-like, whose mass m x o is approximately given by the 
soft Susy breaking gaugino mass ~ M\. Assuming unification of the three gaugino masses 
for the bino (Mi), the winos (M 2 ) and the gluino (M 3 ) at the scale of Grand Unification, 
Mi is naturally the smallest among these mass terms at the electroweak scale. However, 
due to the lower bound of ~ 100 GeV on M 2 from the lower bound on chargino masses, the 
assumption of unification of the three gaugino masses implies Mi > 50 GeV and a similar 
lower bound on the mass of the LSP. 

The assumption of unification of the three gaugino masses can be dropped, however, in 
that case Mi and hence the LSP mass m x o can be arbitrarily small. Then, on the other 
hand it can become difficult to satisfy the WMAP constraint on the dark matter relic 
density i.e., to reduce the dark matter relic density after the Big Bang to an acceptable 
value compatible with this constraint. To this end, dark matter annihilation processes have 
to be sufficiently effective. For a LSP mass > 50 GeV the following pair annihilation 
processes can be relevant: exchange of Susy partners of fermions (sfermions, in particular 
sleptons) in the t-channel, and Z-exchange or Higgs-exchange in the s-channel (if the LSP 
has a sufficiently large higgsino component). In addition, neutralinos can co-annihilate 
with other sparticles if they have similar masses, but co-annihilation processes will not be 
relevant for a light LSP as considered here. In the MSSM, sufficiently effective dark matter 
annihilation processes impose constraints on a light LSP: 

Considering LSP annihilation via slepton exchange in the t-channel, a lower bound 
m x o > 18 GeV was derived in [8, 9] from the lower bound of ~ 100 GeV on the slepton 
masses. (Relaxing this bound to > 80 GeV for stau masses, one obtains m x o > 13 GeV 
[8-10], unless the LSP mass is very small corresponding to hot dark matter [10].) Allow- 
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ing for LSP annihilation via CP-odd Higgs (A) exchange in the s-channel, a lower limit 
m x° ~ ^ GeV was given in [11-14] from m^ > 90 GeV for large values of tan/3 > 25. 
However, as noted in [15], this region of the parameter space of the MSSM is now strongly 
constrained by the bounds on B s — > fi + fi~. A LSP with a mass in the 5 — 15 GeV range in 
the MSSM has been considered in [16] without, however, asking for a correct relic density. 

In the Next-to-Minimal Supersymmetric Standard Model (NMSSM, see [17,18] for recent 
reviews), which can solve the /i-problem of the MSSM, the Higgs and neutralino sectors 
are extended by gauge singlet states. As noticed in [19-22], the mass of the LSP can be 
considerably smaller in the NMSSM than in the MSSM and can still be compatible with 
the WMAP constraint on the relic density. This is a consequence of a light CP-odd Higgs 
boson in the spectrum (on top of the CP-odd Higgs boson of the MSSM), which can be 
mostly singlet-like and which is not ruled out by LEP-constraints. Then, sufficiently large 
LSP annihilation cross sections via the exchange of this additional CP-odd Higgs boson in 
the s-channel may be possible even for a light LSP with mass of a few GeV. 

A light LSP in the NMSSM could be a (dominantly) singlet-like state; in this case, 
however, its direct detection cross sections would be tiny. On the other hand, as in the 
MSSM, a light LSP in the NMSSM can originate from a small value of Mi in which case 
it will be dominantly bino-like and can have larger direct detection cross sections. These 
have been estimated in [21], where also constraints on the corresponding parameter space 
from B-physics, LEP and T-physics were discussed. However, the corresponding points in 
the parameter space given as examples in [21] suffer from a negative effective /^-parameter 
(which is in conflict with the measured anomalous magnetic moment of the muon), and not 
all experimental constraints considered below are taken into account. 

In view of the interest in a light LSP with a mass in the 2 — 20 GeV range, we find it 
appropriate to study upper bounds on its direct detection cross sections in the NMSSM. 
Direct detection cross sections in the NMSSM including WMAP constraints have been 
studied before in [23-29], but not for the LSP mass range considered here. Apart from 
WMAP constraints, we take care of a lengthy list of experimental constraints from B- 
physics (important for large tan (3 and/or relatively light charged and CP-odd Higgs bosons 
as relevant here), T-physics and LEP-constraints on neutralino production. Among the 
latter, OPAL limits on e + e _ — > xlxl (i > 1) turn out to be very important. Since these 
are also relevant for the MSSM, but have hardly been discussed before (a notable exception 
is [10]), we study their consequences in some detail. For the numerical analysis we use 
the code NMSSMTools [30,31] coupled to micrOMEGAs [32-34]. As a result we obtain 
upper bounds on the spin-independent and spin- dependent LSP-nucleon cross sections of 
a SI < 10 -42 cm -2 and a SD <4x 10~ 40 cm -2 , varying somewhat with the LSP mass in the 
2 — 20 GeV range. The maximal value for a SI is indeed near the estimate given in [21]. 

In the next section (2) we present the relevant parameters of the NMSSM and their 
impact on the LSP cross sections. In section (3), we discuss the relevant experimental 
constraints. The consequences of the OPAL constraints on e + e~ — > XiXi (i > 1) on the 
parameter space (implying a lower bound on /i e //) are estimated in an analytic approxima- 
tion, which reproduces well the full numerical results. Section (4) is devoted to our results 
and conclusions. 



2 The NMSSM and the impact of its parameters on 
the LSP cross sections 

In the NMSSM the /i parameter of the MSSM is replaced by a Yukawa coupling A to a 
gauge singlet (super-) field S. Then, the vacuum expectation value (vev) s of the real scalar 
component of S generates an effective //-term 

H eff = XS . (1) 

Most of the time one studies the NMSSM with a scale invariant superpotential W which 
contains, apart from the Yukawa coupling of S to the MSSM-like Higgs doublet fields H u 
and Ha, a trilinear term ~ ^S 3 . Hence the Higgs mass term fiH u Hd in Wmssm is replaced 

by 

Wnmssm = XSH u H d + —S + . . . . (2) 

(Occasionally one considers the so-called nMSSM [22,25,28] without the trilinear coupling 
~ | S 3 , which is replaced by a tadpole-term ~ £,fS.) Compared to the MSSM, the gauge 
singlet superfield S adds additional degrees of freedom to the CP-even and CP-odd Higgs 
sectors as well as to the neutralino sector. Hence the spectrum contains 

• 3 CP-even neutral Higgs bosons Hi, i = 1,2, 3, which mix in general; 

• 2 CP-odd neutral Higgs bosons A\ and A 2 ; 

• one charged Higgs boson H ± ; 

• five neutralinos x°n £ = 1 . . . 5, which are mixtures of the bino, the neutral wino, the 
neutral higgsinos and the singlino; 

• two charginos which are mixtures of the charged winos and the charged higgsinos. 

Apart from the Susy generalisations of the Standard-Model-like gauge and Yukawa 
couplings and the superpotential in Eq. (2), the Lagrangian of the NMSSM contains soft 
Susy breaking terms in the form of gaugino masses M 1; M 2 and M 3 for the bino, the winos 
and the gluino, respectively, mass terms for all scalars (squarks, sleptons, Higgs bosons 
including the singlet S) as well as trilinear scalar self-couplings as \A\SH u Hd + ^A K S 3 , 
which reflect the trilinear couplings among the superfields in the superpotential. 

Expressions for the mass matrices for all Higgs- and neutralino states can be found 
in [17,18]; below we confine ourselves to those which are of relevance subsequently. Dropping 
the Goldstone mode, the 2x2 mass matrix for the CP-odd Higgs bosons M. 2 p in the basis 
(Amssm, Si) has the elements 



2 /ic ff (A x + ks) _ 
sin 2/3 



Mp > u ~ sin 2/3 ~ Ma " 

•Mp,22 = \(A x /s + AK)v u v d - 3kA k s , 

Mp 12 = \{A x -2ks)v (3) 



where v u , Vd denote the vevs of H U1 H d , respectively, v = \Jv\ + v\ ~ 174 GeV and, as 
usual, tan/3 = v u /vd- The matrix element M. 2 Pll would be the mass squared of the MSSM- 
like CP-odd scalar Amssm, if the singlet sector were absent; subsequently we will denote it 
simply by M\. (This parameter can replace the parameter A\.) For any (possibly large) 
value of M\, M. 2 P can have another small eigenvalue corresponding to an additional light 
CP-odd Higgs boson Ai which is mostly singlet-like. This state will be relevant for the LSP 
annihilation cross section below. 

The mass of the charged Higgs scalar is given by 
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M 2 H± =M 2 A + v 2 (^-\ 2 ); (4) 



note that it decreases with increasing A. As is well known, too small values of M#± can 
cause disagreements between measurements and corresponding contributions to B-physics- 
observables as b — > sj; this will be of importance below. 

Notably for large Ma, one of the 3 CP-even Higgs bosons will have a mass close to Ma- 
In the MSSM, the corresponding CP-even state is denoted by H, and we will maintain this 
denomination. The spin-independent LSP-nucleon cross section will be dominated by the 
exchange of this CP-even scalar H, since its couplings to down-type quarks (particularly 
the strange quark) are enhanced for large values of tan/3. 

Also, the mass of the charged Higgs scalar is close to Ma for large Ma', then the states 
H, Amssm and H ± form a nearly degenerate SU(2) doublet. In fact this approximate 
degeneracy holds down to fairly low values of Ma ~ 300 GeV. 

In the neutralino sector, the bino Ai and the neutral wino A 2 mix with the neutral 
higgsinos ipd,ip® and the singlino ips, an d generate a symmetric 5x5 mass matrix .Mo- in 
the basis ip° = (—iXi,—i\l,ip^,ip^,ips), -Mo reads 
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(5) 



It can be diagonalized by an orthogonal real matrix Nij such that the physical masses m x o 
ordered in \m x o\ are real (but not necessarily positive). Denoting the 5 eigenstates by Xi, 
we have 

X° = N^ . (6) 

Finally, the chargino masses are described by a 2 x 2 mass matrix containing M 2 and 
fj, e ff as diagonal entries. The lower bound of ~ 103 GeV on the lightest chargino implies 
at least the constraint 

Mm{M 2 , |/i e// |}> 100 GeV (7) 

(one can choose M 2 > by convention). 

Next, we discuss the dominant contribution to the spin- independent LSP-nucleon cross 
section a SI . Leaving aside scenarios with light squark masses of ~ 100 GeV (which are 



difficult to reconcile with Tevatron constraints), a SI is dominated by the exchange of CP- 
even Higgs bosons, which couple mostly to the strange quark sea. Among the CP-even 
Higgs bosons, the coupling of the state H to down-type quarks (as the strange quark) 
increases with tan/3. Hence, although its mass is generally larger than the mass of the 
Standard-Model-like Higgs boson h, if-exchange provides the leading contribution to a SI 
for large values of tan /?. Then, the dominant component of H is given by H^. 

The dominant coupling of H ~ H^ to the LSP is induced by the bino-higgsino-Higgs 
vertex ~ g\ and hence proportional to giNuN^, where Nu denotes the bino- and A 13 the 
^-higgsino- component of the LSP. All in all one finds 

a** ~ Nl x Nl^ , (8) 



m H 

which shows that the largest values of a SI are obtained for a large product NuNiz, large 
tan/3 and low values of m#. 

The dominant contribution to the spin- dependent LSP-nucleon cross section a SD orig- 
inates, as in the MSSM, from Z-exchange. At first sight one could imagine that, for a 
light CP-odd Higgs boson A l5 its exchange could also give important contributions to a SD . 
However, a light A\ is dominantly singlet-like and, moreover, the coupling of its doublet 
component to strange quarks is always tiny compared to the Z-boson coupling. 

The coupling of the Z-boson to the LSP originates from the gauge couplings of the 
higgsino components ip® and ip®. Since no additional free parameters intervene, the spin- 
dependent cross section a SD is proportional to 

a SD ~ (N? 3 - Nltf . (9) 

Finally the LSP annihilation cross section a ann is dominated, for the LSP mass range 
2 — 20 GeV under consideration, by the exchange of a light A\ in the s-channel. The 
dominant contribution to the A\x\x\ coupling is induced by the doublet component of A 1 
and the bino-higgsino components of Xi as i n the case of the Hxlxl coupling above; the 
singlet components of A\ and Xi pl & y a minor role here. In any case one has (neglecting 
the finite width of A\ and the velocity of Xi near the freeze-out temperature) 



(mi — 4m o 



2 \2 



(10) 



and hence a ann can be sufficiently large for suitable values of m 2 Ai , the lightest eigenvalue 
ofA^inEq. (3). 

3 Experimental constraints on the parameter space 

In this section we discuss various constraints on the parameters of the NMSSM, notably 
(but not exclusively) from LEP and B-physics, separately in various subsections. 



3.1 Constraints from sparticle and Higgs searches 

As we have seen in Eq. (8), a large spin-independent detection cross section a SI requires 
bino- components Nu and higgsino- components Ni 3 of the LSP. For small Mi such that m x o 
is in the 2 — 20 GeV range, the bino component of Xi is automatically large. However, a 
large higgsino component of Xi require relatively small values for fi e ff (below ~ 160 GeV) 
in the mass matrix .Mo in Eq. (5). Consequently the neutralino states x% an d X 3 (f° r 
M2,2ks > fi e ff) are higgsino-like with masses of the order of /i e //- Then, the production 
process e + e _ —¥ XiXi (& = 2, 3) was kinematically possible at LEP2, and corresponding 
limits from DELPHI [35] and OPAL [36] have to be taken into account. 

The strongest limits come from OPAL at 208 GeV, where we can assume 100% Z* 
branching ratios for the x°i decays (see Fig. 10 in [36]). Upper bounds on the cross section 
are given in 5 GeV- wide bins of m v o. Since we will find m v o — m v o ~ 40 GeV, the bounds 
apply for x\ and x% separately. For m x o < 20 GeV, at least one of the X2 or X 3 production 
cross sections (in association with x°) is bounded from above by 0.05 pb. 

In principle, both Z*-exchange in the s-channel and selectron exchange in the t-channel 
contribute to this cross section. However, the interference between these channels is positive, 
hence the most conservative bounds on the parameters are obtained by assuming heavy 
selectrons and that e + e~ — > XiX^ originates from Z*-exchange only. The expression for 
c"z(e + e~ — i x?Xi) is given, e.g., in [37] and can be written as 

a z (e + e- -» X ix" ) = J 9 } + ^L, (N l3 N i3 - N l4 N iA ) 2 (\ - sin 2 6 W + 2 sin 4 6 h 



327r(s-M|) 2V iJ w i4 " 7 \A 

(note the different basis for the neutralinos in [37]) with 

A(s) = s 2 + m\ + m\ - 2s ( mK + m 2 ) - 2m 2 /0 m 2 /0 ■ (12) 

Ai Xi \ Ai Xi J X\ Xi 

In order to obtain an approximate expression for the resulting constraints on the pa- 
rameters, we first neglect m x o everywhere in (11). Using numerical values for the gauge 
couplings, (11) simplifies to 

(s — m 2 () ) 2 / m 2 \ 
a z {e + e~ -► xlx°)W\ ^ 4.9 x 10 4 — _ ^ f 1 + -£ j (N l3 N t3 - N lA N iA f (13) 

with s and the masses in GeV. Next we look for approximations for the relevant neutralino 
mixing parameters Ny. For simplification we assume M 2 , 2ks ^> |/i e //| such that the 
wino- and singlino-sectors in the mass matrix A4 in Eq. (5) decouple. (The wino- and 
singlino- components of the LSP hardly contribute to the spin-independent cross section.) 
Assuming, in addition, large tan/3 such that Vd <^C v u , .Mo can be diagonalised analytically 



with the results (we define u = g\V u j \[2 ~ 43 GeV and write \x = \i e ff) 



N n == , N 13 ;= , N : 



1 + $ */! + *? 



2 V J- t rr 



'11 



iV 2 i ~ — A?3i ~ , /V 2 3 ~ -^33 ~ / , ^"24 ~ ^34 ~ ~ 7= • (14) 

AA + g A/i + £ V2 



Replacing these expressions into (13), using the numerical values for s and Mz in the 
denominator and, notably, approximating mo ~ u, one ends up with 

*z(e + e- -* X?X°)[pb] * 8.3 x 10' 7 ^^ (l + f) > ( 15 ) 

where s, /x and w are in GeV. Then the upper OPAL bound on az{e + e~ — > XiXi) of 0.05 pb 
becomes a lower bound on |/i| (= /x e //), 

|/i e// |>lllGeV. (16) 

A somewhat stronger version of the OPAL bound (er^ < 0.01 pb) is implemented in 
the default version of NMSSMTools [30, 31]. We replace it by the published value of 
0.05 pb [36] for our numerical analysis. From this, without any approximations, we ob- 
tain \fi e ff\ > 114 GeV (varying somewhat with M 2 and tan/3) for small values of m x o in 
good agreement with the previous estimation. We remark that, within the approximations 
used in Eqs. (14), this implies an upper bound on N i3 of Nf 3 < 0.12. 

Next, we consider constraints from the upper bound on the invisible Z decay width, to 
which the decay Z — > x?X? would contribute. From [38] we obtain AT 1 ™ < 2.0 MeV (a 
value slightly above the one used in [10], but below the value used in [21]). The expression 
for the contribution to AT 1 ™ from \\ reads 

ATT = ^ W3 " Nf A ) 2 (l - ^ ~ 0.165 GeV (JV? S - A^) 2 , (17) 

where the last expression holds for small rn x o. Then the upper bound on AT 1 ™ implies 

K 2 3 - N u\ < 0-11 ■ (18) 

For large tan/3, where Nf A <C Nf 3 , this bound on Nf 3 is very similar to the bound obtained 
above from the OPAL limits. According to the numerical analysis without approximations 
we find that the constraints on the parameter space from e + e~ — > X\Xt are mostly somewhat 
stronger than those from AT 1 ™; from (8) and (9) it should be clear, that these constraints 
are relevant for upper bounds on the spin-independent and spin-dependent LSP-nucleon 
cross sections. 

For the chargino masses we require a lower bound of 103 GeV [39], which implies lower 
limits on combinations of the parameters M 2 and fi e ff- I n the neutral Higgs sector we 
apply the various constraints from [40]. Since the lightest CP-even Higgs boson h is mostly 
Standard-Model-like in our case, these constraints reduce to the well-known bound m^ > 
114 GeV. On the other hand the constraints from B-physics, as described below, will imply 
charged Higgs masses above ~ 200 GeV, hence additional bounds from direct charged Higgs 
production are not required. 



3.2 Constraints from B-physics 

Relevant constraints from B-physics originate from bounds on BR(b — > 37) = (3.55±0.24± 
0.09) x 10~ 4 [41], AM, = 17.77 ± 0.12 ps~ x [42] and AM d = 0.507 ± 0.005 ps" 1 [41], and 
the branching ratios BR(B S — > /U + /i~) < 5.8 x 10~ 8 [43] (which was recently improved to 
< 4.3 x 10" 8 at 95% C.L. [44]) and BR(B + -+ t + u t ) < (1.67 ± 0.39) x 10~ 4 [41]. These 
constraints are implemented in NMSSMTools as described in [45], to which we refer for the 
corresponding contributions to these observables in the NMSSM. 

It should be noted that charged Higgs boson exchange contributes to BR{b — > 57) and 
BR(B + —7- r + z/ T ), hence the corresponding limits impose lower bounds on m H ±. On the 
other hand, Susy diagrams also contribute to these observables which depend on parameters 
like M 2 , /i e //, M squark and A top [46]. For specific choices of these parameters (notably not 
too large positive values of A top ), the charged Higgs boson contributions can be partially 
cancelled. This will be relevant below, since the spin-independent LSP-nucleon cross section 
(8) is maximal for small win and, as noted above, m# ~ m#±. 

At large tan/3, the observables AM S , AM^ and BR(B S — >■ /i + /U~) can receive large 
contributions from a light CP-odd Higgs boson Ai [45] which, in turn, plays an impor- 
tant role for the LSP annihilation cross section (10) for a small LSP mass rn x o. Again, 
additional Susy contributions (box diagrams) exist, leading to a complicated combination 
of constraints in the parameter space. We find that, for a small LSP mass (light A\), 
practically all these observables impose bounds on various corners in the parameter space. 

3.3 Additional constraints 

On the dark matter relic density we impose the 3 a WMAP bound [7] 

0.091 < tth 2 < 0.129 , (19) 

which requires a sufficiently large LSP annihilation rate (10). 

A light CP-odd Higgs boson A\ with a mass below ~ 9.3 GeV can appear in radiative 
T —)• ^4x7 decays, on which CLEO [47] and BaBar [48,49] have obtained upper bounds. 
These can be translated into the parameter space (couplings of A\) of the NMSSM [50-52] 
and are implemented, together with constraints from possible A\ —rjb mixing effects [51], in 
NMSSMTools. We find that these constraints are so strong (imposing, essentially, strong 
upper bounds on the A\bb coupling for rriA 1 ^10 GeV) that it becomes very difficult to 
obtain a LSP annihilation rate compatible with (19) for m x o < 2 GeV. 

Finally we require that the Susy contributions to the anomalous magnetic moment of 
the muon (see [53, 54] for such contributions in the NMSSM) improve the disagreement 
between the result of the E821 experiment [55] and the Standard Model; as in the MSSM, 
this implies a positive value for fi e ff- 

4 Results and conclusions 

Before we turn to our results, we discuss the range of parameters used to maximise the 
direct detection cross sections respecting the experimental constraints above. First, for A 
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we choose a small value A = 0.05 such that its negative effect on M H ± as in Eq. (4) remains 
negligible while a non-zero doublet component of A\ induced by the off-diagonal term in 
Eq. (3). A large value for k = 0.55 makes the singlino (and the singlet-like CP-even Higgs 
state) heavy such that perturbing mixing effects in these sectors are avoided. 

For the Susy breaking squark and slepton masses we use 1 TeV such that sleptons 
hardly contribute to e + e~ — > XiX?- Map varies from 300 to 650 GeV where if ± -induced 
and Susy-induced contributions to BRib — \ 37) tend to cancel [45]. The Susy breaking 
gluino and the wino masses are chosen as M 3 = 350 GeV and M 2 = 180 GeV, respectively. 
(These parameters appear in the loop-induced flavour changing Ax-quark vertices [46], 
which should be small in order to allow for a light A\ consistent with the constraints from 
BR(B S -»• fi+fi-).) 

Although Eq. (8) suggests that a SI is maximised for very large values of tan/3, the best 
compromise with B-physics is obtained for reasonable values of tan/3 ~ 35 — 44. Likewise, 
Eq. (14) suggests that /i e // should be as small as possible in order to maximise NnNi 3 , 
but we find that the best compromise in parameter space is obtained for a e ff ~ 128 GeV. 
Eq. (8) also suggests that a SI is maximised for rrin as small as possible. However, we 
recall that win ~ M4 ~ rnH± and that m#± is bounded from below by several B-physics 
processes. We choose Ma as an input parameter of the NMSSM (instead of A\) and find 
the largest values of a SI for Ma ~ 260 — 315 GeV, implying ran ~ 205 — 260 GeV and 
mn± ~ 225 — 280 GeV where the larger values correspond to lower masses of m x o below. 

The Susy breaking bino mass term Mi determines m x o, with Mi ~ 23.5 GeV for m x o ~ 
20 GeV and M x ~ 3.0 GeV for m x o ~ 2.0 GeV. Finally A K is chosen in the range A K ~ 
—14 ... — 4 GeV, which determines rriA 1 such that Xi P arr annihilation via A\ exchange 
in the s-channel gives the correct relic density in agreement with the WMAP bound in 
Eq. (19). Due to the relatively large couplings involved in the Xi P arr annihilation process, 
A\ must actually be off-shell and hence rriA l substantially larger than 2m x o; otherwise the 
relic density is too small. In Fig. 1 we show the result for rriA l as function of m x o . 

For itia 1 < 40 GeV ijn x a < 5 GeV) the constraints from BR(B S — > fi + fi~) (where 
Ai appears in the s-channel) become particularly strong and require a somewhat smaller 
doublet component of A\. Denoting its doublet component by sin^, we have sin 9a ~ 0.8 
for m Al > 50 GeV, but sin6U ~ 0.45 for m Al ~ 10 GeV. We note that for m Al < 40 GeV 
the value of A K has to be chosen within a precision less than 1% such that the relic density 
of Xi is below the WMAP bound (possibly smaller), but m 2 A > 0; hence this region in 
the parameter space requires considerable fine tuning. For rriA 1 < 10 GeV (m x o < 2 GeV) 
the constraints from CLEO and BaBar become so strong that sin 9a must be much smaller 
requiring an even stronger fine tuning of parameters, therefore we will not consider this 
range of parameters subsequently. 

The components of Xi (the coefficients Nu, see Eq. (6)) hardly change in the range 
m x o = 2 — 20 GeV considered here, once we maximise the product N11N13 in order to 
maximise a sl . We have 

N u 0.94 , N l2 ~ 0.01 . . . 0.03 , /Vi 3 ~ -0.32. . . - 0.34 , 

N14 ~ 0.013 . . . 0.06 , V 15 ~ 0.001 . (20) 

The masses of the mostly higgsino-like states X2 an d X3 are ~ 105 and ~ 145 GeV, respec- 
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Figure 1 : vfiA x as function of m x o such that the relic density of Xi is in agreement with the 
WMAP bound Eq. (19). 



tively, and hence as stated before, the limits on az{e + e~ — > XiXl) are relevant. 

The scattering rates of Xi depend somewhat on astrophysical parameters as the escape 
velocity v max and the dark matter density po near the sun and, more importantly, on nuclear 
form factors (quark matrix elements) as the pion-nucleon sigma term u^n and the size of 
SU(3) symmetry breaking parametrized by cr . (The difference a — a^N is proportional 
to the strange quark matrix element.) For the astrophysical parameters we use the default 
values of micrOMEGAs v max = 544 km/s and po = 0.3 GeV/cm 3 [34]. The default values 
in micrOMEGAs for a^N and Cq are u^n = 55 MeV and er = 35 MeV. 

The corresponding results for the upper bound on the spin-independent cross section of 
Xi off protons a^ 1 in the NMSSM are shown in Fig. 2 as a function of m x o as a full red line. 
(The spin-independent cross section off neutrons is nearly the same.) In order to indicate 
the variation of this upper bound with a n N and o"o, we show a red dashed line as the upper 
bound on a^ 1 for a n N — 73 MeV and a = 30 MeV, which would correspond to a larger 
strange quark matrix element and hence an increase of a SI by a factor ~ 3.3. 

Also shown in Fig. 2 are the regions compatible with the excesses of events reported by 
DAM A [1] (without channeling (dark blue) and with channeling (light blue)), CoGeNT [2] 
(light green) and a fit to the two events observed by CDMS-II [56] (denoted as CDMS- 
09 fit surrounded in dashed green; these events are also compatible with background). 
Exclusion limits are shown from XenonlO [4] (violet), XenonlOO [5] (black) and CDMS- 
II [3,6] (magenta, assuming that the two observed events originate from background). 

Fig. 2 is our main result, which leads to the following conslusions: 

• It seems difficult to explain the excesses of events reported by DAMA and CoGeNT 
within the general NMSSM (without unification constraints on Mi). Hence, as stated 
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Figure 2: Upper bounds on the spin-independent cross section a^ 1 in the NMSSM for 
default values of the strange quark content of nucleons as a full red line, and an enhanced 
strange quark content of nucleons as a dashed red line. Also shown are regions compatible 
with DAMA, CoGeNT and CDMS-II, and limits from XenonlO, XenonlOO and CDMS-II 
as explained in the text. 



in [21], significant modifications of parameters like a larger local dark matter density 
Po would be required to this end. On the other hand, the two events observed by 
CDMS-II (within the contour denoted as CDMS-09 fit) could be explained in the 

NMSSM. 

• Actual limits of XenonlO, XenonlOO and CDMS-II on spin-independent cross sections 
of WIMPS in the 2 — 20 GeV mass range test regions of the parameter space of the 

NMSSM. 

For completeness we have also considered the spin- dependent cross section a SD in the 
NMSSM, which is maximal for tan/3 > 20 (such that iVf 4 < N? 3 in Eq. (9)), large values 
of Ma (since m# is irrelevant here), and /i e // ~ 121 — 129 GeV. In Fig. 3 we show the 
maximum of the spin- dependent cross section off protons aZ for the same range of m x o = 
2 — 20 GeV. Note that a SD originates from Z-exchange, hence the spin-dependent cross 
section off neutrons a^ D is given by a^ D ~ 0.78 x a^ . The actual experimental upper 
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limits on a SD are one to two orders of magnitude larger [57] than the upper bounds in the 
NMSSM and not shown in Fig. 3. 
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Figure 3: Upper bounds on the spin-dependent cross section or in the NMSSM 



To conclude, we have performed a detailed analysis of the parameter space of the 
NMSSM for general values of M l5 which allows for WIMP masses in the 2-20 GeV 
range. In contrast to the MSSM, light bino-like WIMPs can have a relic density compatible 
with WMAP constraints due to a light NMSSM-specific CP-odd Higgs state which can be 
exchanged in the s-channel. Due to reported excesses of events compatible with WIMP 
masses below 20 GeV, this region is of particular interest. 

We have studied in detail the constraints on this region of the parameter space of the 
NMSSM from LEP and B-physics, and the regions of parameter space which give rise to 
maximal direct detection cross sections while not contradicting experimental limits. The 
resulting upper bounds on <j SI < 10~ 42 cm 2 = 10 -6 pb make it difficult to explain the 
excesses of events reported by DAMA and CoGeNT within the NMSSM for small values 
of M\. On the other hand, the two events observed by CDMS-II could be explained in the 
NMSSM. 

Notably the XenonlO limits [4] on a SI for WIMP masses below 20 GeV start to test 
corresponding regions of the NMSSM parameter space. Future results from XenonlOO could 
confirm the presence of a light WIMP compatible with the NMSSM, or impose further 
constraints on its parameter space. 
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